Operator-sum representation of time-dependent density operators and its applications 
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We show that any arbitrary time-dependent density operator of an open system can always be 
described in terms of an operator-sum representation (Kraus representation) regardless of its initial 
condition and the path of its evolution in the state space, and we provide a general expression 
of Kraus operators for arbitrary time-dependent density operator of an TV-dimensional system. 
Moreover, applications of our result are illustrated through several examples. 
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Arbitrary states of any quantum system, be it open or 
closed, can always be described by density matrices p{t). 
If the system is closed, its time evolution is unitary, and 
there exists a unitary operator U(t), such that p(t) — 
U (t) p(0)U (t)t , where p(0) is the initial state. However, 
if the system is open, the evolution is not necessarily 
unitary and the above relation between p(t) and p(0) is 
in general not valid. To describe the evolution of an open 
system, one usually employs the Kraus representation^]. 
If there exist operators M M (t), called Kraus operators, 
that satisfy 

p(t)=J2M»(t)p(0)M li (tf, (1) 

and 

J2M^M^t)=I, (2) 

then the evolution of p(t) is said to have operator-sum 
representation or the Kraus representation. 

It has been known that a completely positive map al- 
ways possesses the operator-sum representation while a 
general positive mapmay not do[J H H U H Q, H ©, 
UJ. The papers [ljj, [12], [13( furthermore pursued the ex- 
istence of the Kraus representation by investigating the 
role of initial correlations between the open system and 
its environment, and showed that a map based on the 
reduced dynamics, in general, cannot be described as the 
Kraus representation in the presence of the initial correla- 
tions because an additional inhomogeneous part appears. 
On the other hand, we know that for a qubit system, 
although a non-completely positive map $ : p(0) —* p{t) 
does not possess the Kraus representation with the Kraus 
operators independent of p(0), the state p{t) = $(p(0)) 
can still be expressed in the form of Eg. ffl) with the 
operators M^t) selected for p(0)@, Q, Q- Arbitrary 
density matrices p{t) of a qubit system can always be 
described in terms of the operator-sum representation. 
This start us to wonder whether the same property is 



valid for A-dimensional systems. That is, for any ar- 
bitrary time-dependent state p(t) and initial state p(0) 
of an TV-dimensional system, is it true that the opera- 
tors M^(t) can always be found so that and @ are 
fulfilled? This is a interesting issue as it extends the no- 
tion of the Kraus representation associated with a map 
to that associated with an evolution of the state and we 
find that it is useful to some physical problem. 

In this paper, we prove that the time-dependent den- 
sity operator of an arbitrary A-dimensional quantum sys- 
tem can always be described in terms of the Kraus repre- 
sentation and provide a general expression of Kraus op- 
erators for the A-dimensional system. Some applications 
of our result are also discussed. 

Firstly, we prove that any arbitrary two density ma- 
trices pa and pb of an A-dimensional open system can 
always be connected by Kraus operators. 

Suppose 

JV 

p A =J2pf\A i )(A i \ = U AP d A Ul (3) 

i=l 
N 

PB = Y,Pf\ B i)( B i\ = UbPbUI (4) 

where pf [pf) and \Ai) (\Bi)) are the eigenvalues and the 
orthonormal eigenvectors of the density matrix pa (pb) 
respectively, p d A (p B ) is the diagonal matrix with the en- 
tries pf (pf), and U a (Ub) is the unitary transformation 
matrix, the i-th column of which is just the vector \Ai) 
(\B.i)). We want to prove that there exist the operators 
Mfj,, satisfying 

Pb =J2 M »PaMI (5) 
£-MtM M = 7. (6) 

To find the required Kraus operators, one may write 
as TV x A matrices with unknown elements and then 
directly solve (J5J and © to determine the matrices. 
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However, this method is not feasible for high dimensional 
systems due to the computational complexity. To over- 
come the problem, we first seek diagonal matrices, which 
are generally easier to manipulate. To this end, we look 
for operators M[ such that 



E 



(7) 
(8) 



To solve the above equations, we obtain the matrix op- 
erators M' (fi = 0, 1, 2, ...N — 1) with the entries 



■IJ,+N 



(9) 



where i, j = 1,2, ...iV[15|. 

With the matrices M' we can construct the Kraus 
operators M M satisfying JSJ) and ©. To see this, we use 
the relations, 



p d A = U\ PA U A , 
p% = uy B U B . 
Substituting (JHIJ and {TJJ into (0), we have 
ulp B U B =J2K U aPaUaM'J, 



that is, 



Let 



PB = J2 U B MlU A p A {U B M'^J. 



M„ = U B M'U 



t 

A' 



(10) 

(11) 

(12) 

(13) 
(14) 



it is obvious that defined by lfH|) satisfy JSJ) and JJJ. 

They are the Kraus operators connecting pa with p B . 
Since Ua, Ub and can be deduced directly from pa 
and pb-, the explicit expression of M M is obtained for 
arbitrary density matrices. 

Secondly, having proved that any arbitrary two density 
matrices pa and ps of an TV-dimensional open system can 
always be connected by Kraus operators, we now revert 
to our discussion concerning the time evolution of open 
systems. We show that the time-dependent state of an 
arbitrary open system can always be described in terms 
of the Kraus representation. 

We replace pa and ps in the above demonstration by 
p(Q) and p(t) respectively. Note that for any evolution 
of an ./V-dimensional open system, the state p(t) and the 
initial state p(0) can always be expressed as 



JV 



P(t) = Y,Pi{t)\lpi{t))(lpi{t)l 
i=l 
N 

rto) = £ft(°)M°)Xifc(o)|. 



(15) 
(16) 



By comparing (JTJJl l(TB|) with J3J @ and using JT3J, one 
can immediately write down the Kraus operators M^(t) 



M^t) = U(t)M'(t)U(0)\ 



(17) 



where U(t) (U(0)) is the unitary transformation matrix 
that diagonalizes p(t) (p(0)) given explicitly by 



u(t) =(in{t) Mt) 
u(o) =(Vi(o) Mo) 



M(t) ) , (18) 
iMO)), (19) 



and M' At) (p = 0, 1, ..JV - 1) are given by 



Mm 



( VpiW 

y/pW) 



V o 







M[(t) = 











o Vm(t)/ 



o Vp^tj 



... y/p N -l(t) 





( 

V?^t) o o 
o Vp^) o 





(2 
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The operators M^t) defined by pT jl -pU )! satisfy |TJ) and 
(J2J), giving the Kraus representation of p(t). Generally 
speaking, one cannot assert that any arbitrary time- 
dependent density operator possesses the Kraus repre- 
sentation simply from the result that any two states can 
be connected by Kraus operators. However, here we can 
make the assertion because Pi(t) are always positive and 
(|2"UI is always valid independent of time t. 

So far, we have shown that the time-dependent state of 
an open system can always be described in terms of the 
Kraus representation. A general expression of Kraus op- 
erators for any arbitrary iV-dimensional system is given. 
When the evolution of the open system is given by p(t) 
and p(0), the Kraus operators can be written immedi- 
ately as M M (t) = C/(i)A/;(i)t/(0) t , where U(t), U{0) 
and M'^(t) are explicitly given by {TBI, (JTHJ and J201- 
Since the Kraus operators are not unique, the expres- 
sion l|17fl is only one set of them. The other equivalent 
expressions of the Kraus operators can be obtained by 
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M^it) = ^M^if^Vfxv, where V^„ are the elements of an 

arbitrary unitary matrix V. 

Thirdly, in the following paragraphs, we will illustrate 
some applications of our result. 

The result can help to clarify some ambiguous con- 
cepts. For example, the Kraus representation theorem, 
recalling the well-known representation theorem, which 
states that a map has the Kraus representation if and 
only if it is linear, completely positive (CP) and trace 
preserving[l|, one may thought that the expression p' = 
$(/?) has no Kraus representation if the map $ is not CP. 
However, our result shows that any two states p and p' 
can always be connected by the Kraus operators irrespec- 
tive of the form of the map acting on p. This shows that 
although non-CP map does not possess the Kraus repre- 
sentation, the expression p' = $(p) can still be cast into 
the Kraus representation. For instance, the transposi- 
tion operator T : p — > p T is positive but not completely 
positive. According to the Kraus representation theorem, 
there does not exist the Kraus representation of the map 
T. However, the expression p T = T(p) can still be de- 
scribed by Kraus operators. Let pa = p and ps = p T , 
the Kraus operators can be obtained directly using (|14|l . 
The important key for clarifying the ambiguity is to note 
that there is a difference between the Kraus representa- 
tion of a map and the Kraus representation for a state 
under the action of a map. The Kraus operators describ- 
ing the Kraus representation of a map are independent 
of the state p while the Kraus operators describing the 
Kraus representation of a state p under the action of a 
map may be dependent on the state. A non-CP map has 
no Kraus representation, but p' — %(p), be $ a CP map 
or not, always has the Kraus representation. 

An important corollary of our result is that there al- 
ways exists a CP map between any two quantum states 
and the map can be represented through N Kraus oper- 
ators. This corollary is obvious, because any two states 
can be connected by N Kraus operators M M and the CP 
map can then be defined by the N Kraus operators. The 
corollary shows that even for such two states p and p' , 
where p' is obtained from p by a non-CP map $ : p — > p' , 
one can still find an alternative map $ that is completely 
positive, satisfying $ : p — > p' . 

The conclusion that an arbitrary time-dependent state 
p(t) can always be described by the Kraus representa- 
tion could have other deep applications in the study of 
open systems, especially in the field of quantum infor- 
mation. For example, it is useful for the study on ge- 
ometric phase. As we know, geometric phases of both 
pure state and mixed state under unitary evolutions have 
been clarified[ll, [13, EE El, EH - A new issue is on 
the geometric phases for open svstems under nonunitary 
evolutions. Some papers ja Lil l2_| just use the Kraus op- 
erators M M to define and calculate the geometric phases 
of open systems. As described in ref. |23(, the relative 



phases are defined by a M = argfr[M /1 (r)/9(Q)], and the ge- 
ometric phases can be calculated by making polar decom- 
position of M^{t), such that M^(t) — h^tyu^t), where 
hfity) are Hermitian and positive, and u^t) are unitary. 
The relative phases a M will then lead to the geometric 
phase when u^(t) satisfy the N parallel transport con- 
ditions (Vii(0)|u M (t)+u M (t)|Vi(0)> = 0, (t = 1,2, ...,7V). 
With the scheme showed in the present paper, any time- 
dependent density matrix can be easily written as the 
Kraus representation. So one can transplant the notion 
of geometric phases defined by using Kraus operators to 
arbitrary evolutions of quantum systems by writing p(t) 
as the Kraus representation. 

Another example of its applications concerns the in- 
verse problem of the evolutions of open systems. Sup- 
pose that the density matrix of a system is given as 
a time-dependent function p(t), one wants inversely to 
deduce the evolutional operators or Hamiltonians that 
evolve the initial state p(0) to the state p(t). This is 
an important issue because physicists sometimes need 
to prepare experimentally a quantum system that is ex- 
pected to evolve along a given path in the projected 
Hilbert space. The approach for solving the problem 
is to constitute a closed system by combining the open 
system with an ancilla. The open system will undergo 
nonunitary evolution while the combined system evolves 
unitarily as g(t) = U sa (t)g(0)U sa (t) + , where U sa {t) are 
unitary operators acting on the combined system and 
g(0) = p(0) <g> |0 Q )(0 a | is the initial state. The issue be- 
comes to finding the unitary operators U sa (t) so that 
tr e g(t) — p(t). It is quite a difficult problem in gen- 
eral. However, the present paper can provide an effec- 
tive approach to obtain the unitary operators. As con- 
cluded above, evolutions of an open system always have 
the Kraus representation and the Kraus operators can 
be directly deduced by the given p(t). With the Kraus 
operators M„(t), one can easily obtain the unitary op- 
erator U sa {t). In fact, in order to satisfy tr e g(t) = 
p(t), the elements of U sa (t) in the bases |^i(0)) <g> \j a ) 
are only required to be [U sa (t)}ij,ko — [Mj{t)]ik while 
[U sa {t)]ij,ki (I ^ 0) are arbitrary but keeping U sa {i) to be 
unitary and J7 ao (t)| t =o = I, where \j a ) (j = 0, 1, ...K - 1) 
are the bases of the ancilla. Obviously, there are infinitely 
many such unitary operators U sa {t), so do the Hamilto- 
nians H = iU(t)U(t) + . One can choose the suitable ones 
which can be easily performed in the laboratory. 

It is also interesting to note that the Kraus operators 
provided in the present paper have some intriguing prop- 
erties. Since these Kraus operators are dependent on the 
eigenvectors of "input state" pa but independent of its 
eigenvalues, all the input states with distinct eigenvalues 
but the same eigenvectors will yield the same "output" 
state pb- That is, the map defined by the Kraus op- 
erators can transform all the input states whose Bloch 
vectors lie on the same diameter of the Bloch sphere to 
the same output state. Specially, one may let pa = jjl, 
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and then has Pb — jfJ2 M^M^. We suppose that these 
properties of the Kraus operators may be useful in the 
study of open systems. 
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